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The well-known treatment of asymptotically flat vacuum fields is adapted to pure radiation fields. 
In this approach we find a natural normalization of the radiation null vector. The energy balance 
at null infinity shows that the mass loss results from a linear superposition of the pure and the 
gravitational radiation parts. By transformation to Bondi-Sachs coordinates the Kinnersley photon 
rocket is found to be the only axisymmetric Robinson- Trautman pure radiation solution without 
gravitational radiation. 
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> ■ I. INTRODUCTION 

o ■ 

, Pure radiation or null dust fields are fields of massless radiation which is considered as the incoherent superposition 
, of waves with random phases and different polarizations but with the same propagation direction. The radiation can 
I [ arise from fields of different types, from electro-magnetic null fields, massless scalar fields, neutrino fields or from the 
high frequency limit of gravitational waves. The energy-momentum tensor of pure radiation is 



Tab = Vn-anb, narf^ = 0, ?7 > 0, (1) 

where n° is the tangent vector of the null congruence the radiation propagates along and ry denotes the radiation 
density. The form of is the result of an averaging process, i. e. the equations of the originating fields (Maxwell 
^ equations or Weyl equation) need not to be satisfied. However, (|^) can be derived from a variational principle 
qh, The local conservation law T'^''.^ — implies that the radiation propagates along geodesies 

bi)! n^-.bn^ oc n". (2) 



The energy density rj is not uniquely defined by the decomposition (|^) of the energy-momentum tensor, is un- 
changed under the following rcnormalization of and rj 

fi^^fn" and fl = vr^- (3) 

This annoying arbitraryness in the choice of n° even can't be avoided by an affine parametrization of n°, n°'-b'n^ — 0, 
the vector still might be multiplied e. g. by a constant factor. 

Asymptotically flat pure radiation fields can describe the exterior of isolated radiating sources. The prototype of 
an asymptotically flat pure radiation field, the Vaidya solution is often taken as the exterior field of spherical 
symmetric fluid distributions. More general pure radiation fields, e. g. the pure radiation solutions of Robinson- 
Trautman type, correspondingly could be models for the exterior fields of more general isolated sources. They are 
expected to contain not only pure but also gravitational radiation. 

Gravitational radiation in asymptotically flat space-times can be treated in an elegant way by using Bondi-Sachs 
(BS) coordinates, and the field equations can be solved by series expansion in terms of the luminosity parameter 
[Q. In the next section we will apply this method to the general axisymmetric asymptotically flat pure radiation field. 
In order to solve the field equations step by step in terms of power series as in the vacuum case, the propagation 
vector n°- of the null radiation is adapted to the luminosity distance by fixing n° at We calculate the leading 
order of the radiation density rj and the total amount of energy which is radiated away by the null fiuid. The energy 
balance at J'^ is studied. 
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In the third section the general results are applied to axisymmetric Robinson- Trautman pure radiation solutions. 
By means of a transformation to BS-coordinates in terms of power series we can fit the pure radiation field to the 
luminosity distance. The general expressions for essential quantities, Bondi mass, news function and the leading term 
of the pure radiation energy density are given and we will show that the photon rocket of Kinnersley is the only 
solution of this class which does not contain gravitational radiation. 

Notation: All quantities related to BS-metrics are denoted by capital letters, e. g. coordinates X"- and metric 
components Gab- The metrics have signature (+, +, +, — ). 



II. ASYMPTOTICALLY FLAT PURE RADIATION FIELDS 

BS-coordinates [Q,^, R,U] can be introduced in the region far from isolated sources near J'^ in asymptotically flat 
space-times. For that radially light rays are considered. The corresponding null congruence is geodesic and affinely 
parametrized, expanding and non-rotating, such that the tangent vector is a gradient, 

W"=0, i^bt^'^O, e.a^O ta = -Ua- (4) 

The function U is taken as null coordinate and the intersection points of the null geodesies with a null hypersurface 
f/=const. are parametrized by the spherical coordinates < < tt and < $ < 2tt. The distance between such an 
intersection point and the source is measured by the luminosity distance R which is normalized such that the closed 
two-surfaces i?=const., C/=const. have the surface of a sphere with the radius R. Near , that is for sufficiently 
large R, R > Rq, the line-element in BS-coordinates has the form 

ds^ = R^hAsidX^ - W^dU){dX^ - W^dU) - 2e^f'dUdR + -^e^^dC/^, (5) 

it 



with 



2hABdX^dX^ = {e^'' + e^^)de^ 4 sin 6 sinh(7 - (5)ded$ + sin^ 9(6"^'^ + e"2*)d$^ (6) 



Here P,j,5,V und W^, A,B e {1,2}, are functions of the coordinates O, R and U , because of the axial symmetry 
they are independent of <i>, Gafc,$ — 0. The asymptotical behaviour of the metric (^) is given by the boundary 
conditions. The metric functions can be expanded in powers of 1/R in the coordinate range O<0<7r, 0<$< 27r, 
Ui <U <U2 and R > i?o, 

lim — = -1 and lim RW^ ^ lim /3 = lim 7= lim (5 = (7) 

R — >oc R R — >oo R — *oo R — *oo R — >oc 

hold. The line element (||) and the boundary conditions (|^) are conserved under transformations of the Bondi-Metzner- 
Sachs (BMS) group §, g. 

We are going to consider space-times containing outgoing pure radiation. In general, the geodesic null congruence 
the pure radiation is propagated along is not the null congruence (^) of the coordinate system. In particular, the 
radiation geodesies may have twist. But the asymptotic behaviour of these null congruences is the same, both end up 
at (see figure |l|). The null vector has only a radial component, i° — > (0, 0, 1, 0) for i? — > 00. The tangent vector 
of the radiation null geodesies is also dominated by its radial component, we choose the normalization factor such 
that n" (0, 0, 1, 0) for i? ^ c» holds too, 

n'^ = t^ + o(^^y (8) 

With this normalization the propagation of the pure radiation is fitted to the luminosity distance and the leading 
term of the radiation density 77 is fixed. 
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FIG. 1. Asymptotically flat fleld with pure radiation 



Now we will solve the field equations in terms of power series of 1/i? as it has been done in the vacuum case. For 
relevant results the ansatz 
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with 



i? i?2 ' i?3 ' \^^4 

n21 7i22 „23 / 1 



3 ^ n-^^ n''^2 jj33 / 1 



(9) 
(10) 

(11) 



V ^-R + 2M+^+0[^] (12) 

1 

1 



W-'-^ + ^ + ^ + ^l^) (13) 
^^ = ^ + ^ + ^ + 0(4), (14) 



A3 = 2(5 + 7i,e<^i + (5i,e7i + ^(7i,e7i + '^i,e'5i) + 4cot e7i(5i (15) 
B3 = 25 + 7i,e7i - ^i.e-^i + cot 6(72 - 6l). (16) 

for the metric functions is sufhcient. The regularity at the axis demands 7, S, 7^9 E^nd 5^q to vanish for G = and 
= vr. Additionally, the ansatz 

— + ^ + ^ + (17) 

n2 = 



' i?2 ' ■^h^) (18) 



for the contravariant components of the propagation vector n"" = {71^,71^,71^,71^) is made, which corresponds to the 
normalization (ph. The condition n^n^ = for n° being null is expanded in powers of 1/i? and the coefficients are set 
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equal to 0. This is the way aU equations are solved. We prefer to describe the following steps rather than to present 
the explicit calculations which can easily be reproduced from these expansions. We find 
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(21) 



and the [/-component rt^ of ?i° is fixed 



(ni2)2 



. 2 n 43 

sm B — - — , n 



2 ' - " 2 ' 

One of the ten field equations Rat — KQr/narib determines the radiation density rj, the nine others are to be solved. 
Starting with i?44, we have 



^ = o(— 

i. e. coefficients of powers of 1/R lower than 2 have to vanish. So from i?33n3^^ = 0{1/R^) we get 

7? + ^? 



/3i = 0, 



P2 
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/33 -0. 



The corresponding conditions for Ri^/nin^ and R23/nin^ yield 



^ _7L^±ii^ _ cot 6(71 and sin 6^2 - -^^^^^rr^ - cot 6(71 - 



(23) 



(24) 



(25) 



The coefficients of 1/R in RnUi ^ and Ri2/Tiin3 vanish for j2,u = = S2.U a-nd with the help of the freedom in the 
BMS-group we can put 
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= 5,. 



From the equation R^i/n^^ — kqT] the leading term of the radiation density is calculated 



-0V = ^ 



-2M^u ~ ll.ij - 52 



1 



2 sine 



1 



sine 



(sin2e[7i + 



(26) 



(27) 



For the remaining nine quantities Rab/nanjj (no summation) the l/i?^-term is set equal to kq??. Only eight of these 
nine equations are independent (the ninth equation corresponds to the null condition for the function n"*^), they give 
us a system of algebraical equations for the functions (7, A^, B4, P4 and the [/-derivatives Q,(7, S^u, 73,(7 and S^^ij. 
This result has to be interpreted in analogy to the vacuum case. The system is determined (up to the treated order 
of 1/R) by the functions M, S, Q, 73 and ^3, given on a hypersurface U^const., and by the functions 71 and Si 
and the radiation density rj, which must be given for all U. Note that the temporal development of the mass aspect 
M is determined not only by the news functions 7i_(7 and Si^u but also by the energy density rj of the emitted pure 
radiation. 

This fact becomes even clearer if the energy balance of the system at J'~^ is taken into consideration. The total 
amount £ of emitted pure radiation energy is calculated by integrating the density 77 over the closed surface i?=const., 
C/=const. in the limit R — > 00. With the surface element R^ sinBdOd^ + 0{R) it follows 



£ = lim 2n / r7i?^sin0de. 



(28) 



Introducing the notations (/) for the average of / over the intervall < O < vr and Ai — (M) for the Bondi mass, we 
get from ( p7| ) with the regularity properties of 71 and Si and the convention kq — Stt 

1 



M 



(29) 
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The mass loss results from a linear superposition of pure and gravitational radiation, provided that the propagation 
vector of the pure radiation is normalized such that holds. The corresponding result follows for the linear 
momentum. So for asymptotically flat pure radiation fields the decomposition of the energy-momentum tensor into 
propagation vector n"" and energy density rj is not arbitrary, the normalization is the natural relation between the 
physical and the geometrical fields n° and 



III. ROBINSON-TRAUTMAN PURE RADIATION SOLUTIONS 

Robinson- Trautman space-times admit a shearfree, expanding and non-twisting congruence of null geodesies, which 
is a multiple eigen congruence of the Weyl tensor. The field equations for aligned pure radiation fields of this class 
are completely solved (see 1^ §24.3). In the coordinates iy9,r, m] the line element reads 

2 

ds^ = ^{d^^ + sin^ Mip^)- 2dudr + 2Hdu^ , (30) 

where for the axisymmetric case 

P^P{d,u) and H ^r{\nP),^- — + —-^ (31) 

hold. Kjr^ denotes the Gaussian curvature of the surfaces r=const., u=const., which we assume to be distorted 
spheres, 

if = p2 [(InP)^^^ +cot7?(lnP),^ + 1] =: A In P. (32) 

The vector fc° = (0,0,1,0) is the tangent vector of the affinely parametrized null geodesies, if is chosen as 
propagation vector of the pure radiation, Rah — '*o?7^a^fc, the radiation density reads 
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1 P 

-?7i,„ -I- 3m(lnP),„ -t- -AAlnP - — 
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(33) 



If we put P = 1 we arrive at the Vaidya solution for which the coordinates <y9, r, u\ are BS-coordinates and by 
comparison with the components of the BS-metric (||) and the ansatz (^)-(|lJ) we find M. — m and 71 = <5i = 0. 

Now we are going to the general line element ( |30| ) and transform it into BS-coordinates. The procedure is similar 
to that presented in Q, the BS-coordinates are expanded in power series in terms of 1/r. In the axisymmetric 
non-rotating case we set 

■!> = (34) 
By the line element of the two-surfaces r=const., w=const. the ansatz R ^ r/P is suggested, such that we start with 

e = To + ^ + ^ + O V (35) 



^=p+^o + ^ + o(^), (36) 



^ = ^o + ^ + ^+C^(^), (37) 

where all coefficients depend on -d and u. For the determination of these coefficients we step by step make use of the 
properties 

G44^0 = Gl^ G^^^-l + o(-^, G" = o(^-^V G"g22 = (sin^ ep4)-\ (38) 
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the transformed contravariant metric components 



must have. The expansion of G yields 



G^^ = + O (^^ j , (40) 

from which To^„ = follows. Using the freedom in the BMS-group we can set 

To = d (41) 

for simplicity. From 

G34 = -^ + o(^i^ (42) 

we have 

Uo,u = P, (43) 

which fixes the u-dependence of C/q, the remaining freedom in the i?-dependance also corresponds to a transformation 
in the BMS-group. All higher coefficients can now be expressed in terms of Tq = and Uq. The condition G"*^ — 
yields 

Ui = -^PUl^ and U2 = ^P^Ul^Uo^m (44) 

and from G^"* = we get 

Ti = -PUo,^ and T2 = P'^Uo^^Uom- (45) 

Finally the coefficients in the expansion of R are to be calculated from G^^G^^ = (sin^ <dR'^)~^, the leading term of 
that condition confirms the ansatz R ^ r/P and the higher terms give 

Ro = ^[Uo,i)i) + cotmo^^] (46) 

and 

P PIJ 

Ri = ^[P.M - cotiJF^] - -[C/q,^^ + coidUo.-oY + ^^[C/o,^^fl - cotz^L/o,^]. (47) 
4 o Z 

Thus the coefficients in the ansatz (^)-(^) are determined, terms of higher order can be calculated in the same 
manner by a continued expansion of the relevant conditions. 

The ansatz (^5|)-(^^ is sufficient to extract characteristic quantities from the metric components. From (H) and 
( p^ ) we know the component G^^ and its series expansion in terms of which can be converted into a power series 
of 1 /r by means of ( ^6| ) 

.33._|.-..,_,|,,(_L).,_,^i^,,(^). 

Now the mass aspect M can be read off from the 1 / r-expansion of G^^ as a function of d and u 

M = |i + ^ - ^ - ^[25 , + Scot ^5], (49) 
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where B is defined by 



B = 



sind ( P 



P Vsini? 



(50) 



In the non-rotating case ^ — 5 and = hold for the metric functions, such that the news is given by 71. As the 
mass aspect A/, the function 71 is calculated from the transformed metric 
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'34 



'2/3 
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i?4 



-1 



as 



9 / \ 2 

sin I Uoi) 



2r2 



O 



1 



(51) 



(52) 



By differentiation with respect to U, which can be calculated as 

,. dU 

li,u = 7i,u lim 

r^oo (ju 



we get 



2 _ _ sin^ i9 / 

1l,U - — - 4p2 I ^ 



(53) 



(54) 



Finally the propagation of the pure radiation is adapted to the luminosity distance R, i. e. the radiation null vector 
should be normalized such that (0) holds. The transformed null tetrad fc° 



dr 



1 



1 



1 



1 



(55) 



shows that does not fulfil the normalization condition (^). Moreover, the normalized propagation vector is given 
by 



p + o[- 

r 



(56) 



and the leading term of the correspondingly normalized energy density of the pure radiation reads 
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1 
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(57) 



IV. THE PHOTON ROCKET 

Recently the 'photon rocket' solution of Kinnersley has been the subject of articles where it is proven not to 
contain gravitational radiation - ||T^ . With the results of the foregoing sections we are now able (a) to give another 
proof of this fact and (b) to show that the photon rocket is the only axisymmetric Robinson- Trautman solution with 
pure and without gravitational radiation. 

Kinnersley's solution is interpreted as the field of a particle emitting pure radiation anisotropically, and accelerating 
because of the recoil. In the coordinates [1?, ip, r, u] the line element reads 

~ ~ - ( ~ 2m\ 
ds^ = r'^{[M + asinMuf + sin^ ddip^) - 2dudr - [l-2arcosd dw^ (58) 
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where the functions m and a depend on the retarded time u and can be interpreted as the mass and the acceleration 
of the particle, respectively. Introducing a new coordinate d by 



tan- = tan-, A^A{u), A^u = a, (59) 
the line element (^8|) is cast into Robinson- Trautman form with the metric function 

P = cosh A + sinh Acosi?. (60) 



If this function P is put into the formula for the news function ( |54| ) we get 71 (7 = 0, i. e. there is no gravitational 
radiation. Conversely, solving the equation for vanishing news, we get P = b{u) + c{u) cos 1} which can be put into 
the form ( ^0|) by a coordinate transformation which preserves the form of the Robinson- Trautman line element. 

The result that the Kinncrsley rocket is the only pure radiation solution without gravitational radiation is in 
contradiction to pd| ] where the energy balance for a perturbed Kinncrsley metric is calculated. But it confirms the 
hypothesis in Q that more anisotropic pure radiation solutions than the photon rocket would emit gravitational 
radiation. 



V. CONCLUSION 

We have been studying asymptotically flat pure radiation fields in terms of Bondi-Sachs coordinates. In this 
framework we find a natural normalization for the propagation vector of the pure radiation. The energy balance at 
J'^ shows a linear superposition of pure and gravitational radiation. The application to the Robinson- Trautman pure 
radiation solutions confirms that the Bondi-Sachs coordinates are powerful in the treatment of asymptotically flat 
fields, not only for vacuum but also for pure radiation fields. 
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